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Abstract 

Using first-principles calculations, we systematically study the potential energy surfaces and 
dissociation processes of the hydrogen molecule on the Mg(OOOl) surface. It is found that during 
the dissociative adsorption process with the minimum energy barrier, the hydrogen molecule firstly 
orients perpendicular, and then rotates to be parallel to the surface. It is also found that the 
orientation of the hydrogen molecule at the transition state is neither perpendicular nor parallel 
to the surface. Most importantly, we find that the rotation causes a reduction of the calculated 
dissociation energy barrier for the hydrogen molecule. The underlying electronic reasons for the 
rotation of the hydrogen molecule is also discussed in our paper. 
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I. INTRODUCTION 



It is of great scientific importance to study the dissociation of diatomic molecules on 
metal surfaces (dissociative sticking), to meet some intrinsic interests on mechanisms for 
bond breaking and bond formation and origins of activation energy barriersi^^. The dis- 
sociation of the hydrogen molecule (H2) on metal surfaces are furthermore key events for 
lots of technological applications such as hydrogen storage for fuels^*^, hydrogen caused 
embrittlement^, and heterogeneous catalysis^'^. Among all the metals, the Mg(OOOl) surface 
is one of the most studied prototypes, both theoretically and experimentally^'^i^'^'^iii. How- 
ever, some basic problems concerning the dissociation process and energy barrier still have 
discrepancies yet. As early as 1981, N0rskov et al. for the first time studied the dissociation 
of H2 on the Mg(OOOl) surface, observing a molecularly adsorbed precursor state and an 
energy barrier of 0.5 eV for dissociation^. Later, by using the more precise first-principles 
methods. Bird et al. discovered that although the most energetically favorable site for dis- 
sociation is the bridge site as reported in N0rskov's paper, the dissociation energy barrier 
is 0.37 eV, rather than the reported 0.5 eV, and there is no precursor state^. However, the 
discrepancy on the dissociation energy barrier remains ever since then. Vegge et al. reported 
a value of 1.15 eV for the dissociation, after systematically considered the zero point energy 
of H2 within their first-principles calculations^''. Wu et al. reported a value of 1.05 eV by 
employing first-principles calculations and transition state theory-. Meanwhile, Johansson 
et al. recently organized an experiment to study the dissociation energy barrier for H2 at 
Mg(OOOl), and the obtained barrier is 0.6 ~ 0.9 eV^-, which disagrees with all the reported 
theoretical values. So it is still an open question as to the correct dissociation energy barrier 
for H2 at Mg(OOOl). 

On the other hand, it has been recently suggested that steric effects might be important 
during the adsorption and dissociation of diatomic molecules on metal surfaces^- . For the 
adsorption of D2 on Cu(lll), both experimental and theoretical investigations found that 
dissociation of D2 occurs preferentially when the molecule approaches with its bond parallel 
to the surfacei^. Similar dependence on polar angle has also been theoretically predicted for 
the dissociative sticking of O2 on the Al(lll) surface^^. When considering the dissociative 
adsorption of H2 on metal surfaces, steric effects might be more important because of the low 
inertia moment of H2 and correspondingly high possibility to rotate. In fact, it has already 
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been pointed out that during the dissociative adsorption of H2 on the Pd(lll) surface, the 
molecular axis orientation has a drastic effect and low activation barriers are only met over a 
small range of 6 values from parallel the surface^. And for the adsorption on the NiAl(llO) 
surface, H2 molecules rotate abruptly when they are close to the surface, which allows them 
to adopt the orientation that is more convenient for dissociation (i.e., nearly parallel to 
the surface)^. For the adsorption of H2 on the Mg(OOOl) surface, however, this important 
issue of steric effects has not been considered yet. Motivated by this observation, here by 
using the first-principles calculations, we study the rotation of H2 during the dissociative 
adsorption on the Mg(OOOl) surface and its corresponding influences on the dissociation 
energy barrier. We show that the most energetically favorable path for H2 dissociation at 
Mg(OOOl) is fundamentally determined by the steric effect. 

II. CALCULATION METHOD 

Our calculations were performed within DFT using the Vienna ab-initio simulation 
packag The PW91— generalized gradient approximation and the projector-augmented 
wave potentiali^ were employed to describe the exchange- correlation energy and the electron- 
ion interaction, respectively. The cutoff energy for the plane wave expansion was set to 250 
eV. The Mg(OOOl) surface was modeled by a five-atomic-layer slab plus a vacuum region 
of 20 A. A 2 X 2 supercell was adopted in the study of the H2 adsorption since our test 
calculations have showed that it is large enough to avoid the interaction between adjacent 
hydrogen molecules. Integration over the Brillouin zone was done using the Monkhorst-Pack 
schemei^ with 9x9x1 grid points. The calculated lattice constants of bulk Mg (a, c) and 
the bond length of isolated H2 are 3.207 A, 5.145 A and 0.748 A, respectively, in good 
agreement with the experimental values of 3.21 A, 1.62 A and 0.74 A The calculation 
of the potential energy surface was interpolated to 209 points with different bond length 
(dn-H) and height (/ih2) of H2 at each surface site. 

III. RESULTS AND DISCUSSION 

After geometry optimization for the Mg(OOOl) surface, we build our model to study the 
potential energy surface (PES) of H2 on the relaxed Mg surface. As shown in Fig. 1, there 
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are four different high-symmetry sites on the Mg(OOOl) surface, respectively the top, bridge 
(bri), hep and fee hollow sites. After PES calculations, we find that the dissociation barrier 
of H2 at low-symmetry sites is always larger than at high-symmetry sites, proving that high- 
symmetry sites play crucial roles in the adsorption of diatomic molecules, similar to that 
has been observed on the adsorption of oxygen molecules on the Pb(lll) surface^^. And in 
the following, we will only give the results at the four high-symmetry sites. At each surface 
site, an adsorbed H2 has three different principle orientations, respectively along the x (i.e., 
[1120]), y (i.e., [1100]), and z (i.e., [0001]) directions. Herein, we use top-x,?y,2, bri-x, z, 
hcp-x, y, z and fcc-x, y, z respectively to represent the total twelve high-symmetry channels 
for the adsorption of H2 on the Mg surface. 

Throughout our PES calculations, we find no molecular adsorption precursor states for 
H2 at Mg(OOOl), according well with all previous reports, except for the one by N0rskov et 
al. using jellium model^. Our calculated result for the lowest dissociation energy barrier, 
as well as that in other theoretical and experimental reports, is given in Table. I. The 
minimum energy path (MEP) for the dissociation of H2 on the Mg(OOOl) surface is found to 
be along the bri-y channel, which is consistent with all previous first-principles studies. The 
transition state obtained from our PES calculations along the bri-y channel is at the point 
where (iH-H=l-12 A and /ih2=1-16 A, which accords with previous results^i^. However, as 
we will see later. This transition state needs to be modified after considering the rotational 
degree of freedom of H2. 

At the bridge site, however, we find that the total energy of the H2/Mg system is not 
always smaller along the bri-y channel than along other channels. This is a key point in 
this paper. In fact, at large values of H2 height from Mg(OOOl) surface, we find that the 
total energy is smaller along the bri- 2; channel than along the bri-y channel. To show this, 
we plot in Figs. 2(a)-(c) the two-dimensional cuts of the PES along the bri-x, y^ z channels. 
Correspondingly, the minimum energy paths in Figs. 2(a)-(c) are collected and plotted in 
Fig. 2(d). It can be clearly seen that there is a prominent crossing point in the minimum 
energy paths along hii-y and bri- 2; channels, at which the distance of the H2 molecule from the 
surface takes a value of /iH2=l-26 A. Before this crossing, the total energy of the adsorption 
system along the bri- 2; channel is always lower than along the bri-y channel. This finding 
indicates that H2 prefers to orient perpendicular to the Mg(OOOl) surface until it reaches 
the height lower than 1.26 A. After the crossing point, whereas, the system along the bri-y 
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channel has a smaller total energy than along the bri-z channel, and H2 tends to rotate from 
the bri-2; channel to the hh-y channel. 

We then further study the influence of the molecular rotation on the dissociation energy 
barrier of H2. For this we have calculated the total energy by fixing the mass center of H2 
at the bridge site with the height of 1.16 A and the molecular bond length of 1.12 A, while 
allowing H2 to rotate around its mass center in the y-z plane (see the inset in Fig. 3). The 
calculated angle dependence of the total energy is shown in Fig. 3. Clearly, it can be seen 
that the transition state (namely, the saddle point in the PES for H2 dissociation) should 
be the structure where the orientation of H2 is 21° from the x-z plane. This rotation of H2 
results in a 86 meV modification on the dissociation energy barrier. This finding suggests 
that steric factors that has not been considered in previous theoretical calculations might 
be (at least partially) responsible for their discrepancies with experimental measurement. 

Although it has long been explored for steric effects on the dissociative adsorption of 
diatomic molecules on metal surfaces, the specific reasoning has seldom been discussed yet. 
Herein we will try to find the underlying mechanisms on the rotation of H2 during the disso- 
ciative adsorption on the Mg(OOOl) surface, by analyzing carefully the charge distributions 
and electronic interactions along the adsorption process of H2. Figures 4(a) and (b) show the 
difference electron density for the adsorption system with /iH2 to be 2.00 and 1.16 A along 
the hri-y channel, namely, 

Ap = p(H2 + Mg(OOOl)) - p(H2) - p(Mg(0001))], (1) 

where p(H2 + Mg(OOOl)), p(H2) and p(Mg(0001)) are respectively the electron density of the 
adsorption system, the H2 molecule and the clean Mg(OOOl) surface. To calculate Ap, the 
atomic positions in the last two terms in Eq. 1 have been kept at those of the first term. 
Through careful wavefunction analysis, we find that at the beginning of the adsorption 
process, the molecular orbitals of H2 orthogonalize with electronic states of Mg and thus are 
broadened. As shown in Fig. 4(a), the surface electrons of Mg are repelled from the region 
occupied by the H2 bonding electrons due to these orthogonalizations. This interaction has 
also been observed during the interacations of H2 with other metals such as the Al(lll)^^ and 
transition metal surfaces^^. When the H2 molecule is close enough to the Mg(OOOl) surface 
and come to the transition state for its dissociation, electrons transfer from electronic states 
of Mg to the antibonding orbital of H2, which can be clearly seen from Fig. 4(b). We can 
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also see from Fig. 4(b) that the orthogonahzations between electronic states of Mg and the 
bonding orbital of H2 still exist at the transition state. 

In total, the orthogonahzations between molecular orbitals of H2 and electronic states of 
Mg cause repulsive interactions between electrons of H and Mg, and thus will enlarge the 
total energy of the adsorption system, while the electrons transfer from Mg to H2 causes 
attractive interactions between H and Mg atoms and lowers down the total energy. So 
during the adsorption process of H2, the total energy of the system firstly goes up, then 
lowers down. This analysis explains why an energy barrier is needed for the dissociation of 
H2 on the Mg(OOOl) surface. Moreover, both the orthogonahzations and electrons transfer 
are always weaker along the bri-2; channel than along the bri-y channel. Therefore, at the 
beginning of the adsorption process, when no electrons transfer happens, the total energy of 
the system is smaller along the bri-2; channel than along the bri-y channel. And at around 
the transition state, when the H2 molecule is very close to the Mg(OOOl) surface, electrons 
transfer begins to dominate the molecule-metal interaction. So the total energy along the 
bri-y channel is smaller at the transition state. Herein, the rotation of H2 can be seen as the 
result from the different interactions that respectively favors the hri-y and hri-z channels. 
As a result, the minimum energy path for the dissociation of H2 is neither along the bri-y 
nor along the hn-z channels. And the corresponding transition state is the one where H2 
orients 21° away from the x-z plane, as shown in Fig. 3. 



IV. CONCLUSION 



In conclusion, we have systematically studied the PESs for the dissociative adsorption 
of the hydrogen molecule on the Mg(OOOl) surface. Our results accord well with previous 
reports on the direct dissociative adsorption process. More importantly, we have found that 
the hydrogen molecule does not always orient parallel to the surface along the dissociation 
channel with the lowest energy barrier. At large molecular heights, H2 orients perpendicular 
to the surface. When getting closer to the surface, H2 begins to rotate such that at the 
transition state, H2 orients 21° away from the x-z plane, which causes a 86 meV modification 
on the dissociation energy barrier. We have revealed that this molecular rotation is because 
of the two different interactions between H and Mg, i.e., the orthogonahzations between 
molecular orbitals of H2 and electronic states of Mg and electrons transfer from the Mg(OOOl) 



6 



surface to the antibonding orbital of H2. As a final concluding remark, here based on the 
present results, we would like to point out that steric effects are important to understand 
the adsorption behaviors of H2 on metal surfaces. 
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TABLE I: Our calculated minimum energy barrier for the dissociation of H2 molecules on the 
Mg(OOOl) surface, and other reported results. 
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DFT (GGA) 


1.05 


n 


Experiment 
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List of captions 



Fig.l (Color online), (a) The p{2 x 2) surface cell of Mg(OOOl) and four on-surface 

adsorption sites. Here only the outmost two layers of the surface are shown, (b) The sketch 
map showing that the molecule (with vertical or parallel orientation) is initially away from 
the surface with a hight hH^- 

Fig.2 (Color onhne). Contour plots of the two dimensional cuts of the potential energy 
surfaces (PESs) for the H2/Mg(0001) system as a function of the bond lengths (c/h-h) 
and the heights (hii^), with H2 at the bridge site orienting along x (a), y (b) and z (c) 
directions, (d) Minimum energy paths obtained along the three different channels with 
different heights of H2 {hB.2) from the Mg surface. 

Fig.3 (Color onhne). The total energy of the adsorption system with different orienta- 
tions of H2 at the transition state point along the hri-y channel where (iH-H=l-12 A and 
/iH2=l-16 A. The inset depicts the definition of the angle a, in which the grey area and 
blue balls respectively represent the Mg surface and two H atoms. 

Fig. 4 (Color online). The difference electron density for the H2/Mg(0001) system with 
the height of H2 to be 2.00 (a) and 1.16 A (b). Blue and grey balls respectively represent 
hydrogen and Mg atoms. Dark and dashed lines respectively represent plus and negative 
values, i.e., electrons accumulation and depletion. 
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